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ABOVE THE BUCKLING LIMIT* 

By F. Reinit zhuber 

INTRODUCTION 


When thin— walled columns formed from flanged sheet, such 
as used in airplane cons truct ion , are subjected to axial load, 
their behavior at failure varies according t* the slenderness 
ratio , 


On long columns the axis deflects laterally while the 
cross section form is maintained; buckling results. The re- 
spective breaking load in the elastic range is computed by 
Euler's formula and for the plastic range by the Engesser— 
Karman formula. Its magnitude is essentially dependent upon 
the length. On intermediate length columns, especially where 
open sections are concerned, the cross section is distorted 
while the cross section form is pr es erved; t wis t ing failure 
results. The buckling load in twisting is calculated accord- 
ing to Wagner and Kappus (reference l). On short columns the 
straight walls of low— bending resistance that form the column 
are deflected at the same time that the cross section form 
changes — buckling occurs without immediate failure. Then the 
buckling load of the total section computable from the buckling 
loads of the section walls is not the ultimate load; quite 
often, especially on thin— walled sections, it lies considerably 
higher and is secured by tests. Both loads, the buckling and 
the ultimate load are only in a small measure dependent upon 
length. 

The present report is an attempt to theoretically inves- 
tigate the behavior of such short, thin— walled columns above 
the buckling load with the conventional calculating methods. 

*"Beitrag zur Berechnung gedrilckter, dunnwandiger Profile 
oberhalb der Beulgr enze , 11 Luf tfahrtf ors chung , vol. 19, no. 7, 
July 20, 1942, pp . 240-247. 



HACA TM No. 1077 


2 


The cross section of the column selected was a square, which 
by reason of its symmetrical properties lends itself best to 
mathematical treatment. The results are then checked by tests. 


THEORY 


A closed section (hollow cylinder) of elastic material 
(modulus of elasticity E, Poisson's ratio v) the cross 
section of which is a square with sides of length t and 
sheet thickness s, is centrally loaded with a force P in 
longitudinal direction (direction of cylinder axis). On 
placing through the median plane of one of the four identical 
sides a system of coordinates x, y, z (fig. l)» the defor- 
mations u, v, w are in direction of the coordinates t and 
the stresses o x , Oy , and T along the median plane. 

If the mean compressive stress p = P/4sb is lower than 
the critical buckling stress p* the stresses are o x - — p, 

Oy - 0, T — 0 and the deformations u = — x , v = 

+ u~y, w = + v |r tt • Between the shortening e and the 
E M2 

compressive stress p proportionality exists then according 
to the equation p — Be as far as the buckling limit; at 
the latter p* = Ee*, where e* is the critical compression. 

If the mean compressive stress p exceeds the critical 
buckling stress p* , the side walls wrinkle lengthwise, the 
cross section assumes the form shown in figure 1, on which, 
since the right angles remain at the corners, a bulge in the 
one side area corresponds to a dent in the next one and the 
corners remain square. If the short column is of sufficient 
length, the buckles in longitudinal direction have in the 
area not affected by the disturbed end surfaces, the constant 
length l; for great length, l = b, The behavior of the 
profile is therefore periodic with the length l. The prob- 
lem thus reduces to an analysis of one of the four identical 
rectangular sides between two successive junction lines as a 
plate (dimensions bXi). 

According to the theory of the thin elastic plate with 
great deflections (reference 2) the stresses a x , Oy , and t 

in the plane of the plate follow the equations 
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TT dy ° 

= 0 

dx dy 


which are satisfied by the stress function formula 


( la) 


Oy = $ 


yy 


r y ~ ®xx» r - - $ 


xy 


(it) 


the stress function $ and the sag w being related through 
the fundamental equations 


A A $ “ E (w 2 y - v xx Vyy^ 


and 


12(1 — v 2 ) ^ ^ w ~ ®yy w xx ~~ $ xx w yy + 2 ^xy w xy = 0 


( 2 ) 

(3) 


The deformations u, v in the plane of the plate follow 
from 

2 u x = $yy - V$ xx - ^ w x 2 (4) 

E V y = $ xx “ v *yy *“ f w y 2 (5) 

(The subscripts u, v, w, and $ denote the partial de- 
rivatives conformable to the respective variables.) 

The combined action between the four side walls at the 
corners y = ± b/2 stipulates the compliance with the bound- 
ary conditions 

Wy (x, ± - Wy ( X ± l, ± |.) ( 6a ) 

W (*’ ± "l) = v ( x * * t) (6b) 

CT y (^ x * * 2 "^ = ® 


(6c) 
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T (, 

« 

+ b/2 

s f a T i 
-b/2 

Since no exact integration of equations (2) and (3) has been 
published so far, recourse is had to an approximate method by 
Maguerre (reference 2) and the Eitz formula 


, ± ^ = 0 ( 6d ) 

y = — pbs=— ~ (6e) 


ux ny f m\ 

w = T| + f cos — cos \ 7 / 


is chosen for the deflection, with f as the free value 
which closely approaches the form of buckling occurring short- 
ly after exceeding the critical load and satisfies (6a). 
Equation (7) entered at the right-hand side of equation (2) 
gives 

4 4 (°°° + °°° ¥) 
of which the particular integral 



can be posted. Eor compliance with (6b) to (6d) solutions 
periodic in x direction and symmetrical with respect to the 
x and y axis 

= -(sjI + a co S a cos 2« 

+ B Sinh cob (8b) 

are added to the homogeneous equation A A $ - 0 related to 
(2), after which the complete solution for $ reads 
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$ = $(p) + $(*0 (8) 

Equation (8) for $ entered in (11)), (4), and (5) gives 

(O' cos ap. 


Ov = E 
x 8 


_ 4l£ cos 53S. [(A + SB) Cosh £sz. + B 2EZ sinh^] -p ...(9a) 

'y = [» r (t)’ + TV ( A Oosh ^ + ^T 2 B Slnh ^r)] 


2 ttx / nv «. 
cos — s — l 91)) 


T = 


4rr 2 

l 2 


Eu = — u E 


- E 


’ TT 
8 


["(A + B) Sinh + B Cosh J sin ...(9 c) 

t(£> 

— sin — — — sin cos 

2n J 8n l *> / 



~ |^(A (1 + v) + 2B) Cosh QjZ 


+ (1 + v) B 5^2 sinh -ZZ 


„ a_ 2-rrx 

s in — 5 — — P x 


(9ft) 


_L sin 2nz 
2n 5 n b 


Ev ' - v 1 f (t)" 

- E t (t) ( y ~£; Bto ~ 2 S z si ; 008 8in nf) 

+ ^ [(A (1 + v) - B ( 1 — v ) Sinh ^ 

+ (1 + v) B Cosh ^jcos ~~ + vpy (9e) 

The boundary condition (6a) is satisfied by (7); to satisfy 
(6b) , which with ( 7) becomes : , .- • ,. • .... 

v (*• * t) ■ 


(6b« ) 
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and. also (Sc), ( 6d.) , and (6e); tha-t is, of altogether four 
■boundary conditions, only three integration constants A , B , 
and p are available. A suitable choice is therefore neces- 
sary to ensure satisfactory approximate solutions. 


Of the four possibilities 

( 6b ’ ) , (6c), ( 6d) (10a) 

(6b* ) , (6c) , (,6e) (10b) 

(6b» ), ( 6d ) , (6c) (10c) 

(6c), ( 6d ) , ( 6e ) (lOd) 


available for predicting A, B, and p, (10a) and (10b) 
cancel out, because for these the longitudinal borders of the 
plate stipulate simultaneously the conditions CTy = 0 in 

(6c) and v = constant in (6b 1 ), which can be reciprocally 
excluded and therefore do not eliminate the boundary integral 
in equation ( 13 ) . With ( 10c) and (lOd) it is found from (6e) 
by integration of (9a) for a x over the plate width by ob- 
servation of ( 6d) , that the constant of integration p is 
directly equal to the average value of the compressive stress 


P 


P 

4bs 


( 11a) 


The other two unknown constants A and B are then obtained 
conformably to (10c) and (lOd) from (6b), (6d) or (6c) and 
(6d). Prom (6b) and (6d) the result is A = B = 0 •; this case, 
already investigated by Mar guerre on the effective width of 
the stiffened plate in compression (reference 2), is not as 
suitable an approximate solution as those obtainable from 
(6c) and (6d) and which were therefore used for computing A 
and B. 


Equations (9b) and (9c) afford for y = ± b/2 in accord 
with (6c) and (6d) the equations 


A 


/ \S 

( ) . 


Cosh it. -2- + B 

V 


(t ) 2 "i s . inh " t + e f (xf - 

^Sinh u ^ + tt ~ Cosh tt = 0 


A S inh tt + B ( Sinh tt + tt Cosh 
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hence 


A = 


B «* 


+ 



S inh tt~j- + tr-y- Cosh n — ■ 
Sinh 2 tt y + 2 it -k 

Slnh n j 

Sinh2rr-^-+2TT-| 


(lit) 


(11c) 


The constants of integration A, B, and p given hy (ll) 
satisfy (lOd) but not (6b* ) which specifies that the deforma- 
tion v on the longitudinal edge of the plate shall be con- 
stant and eq.ua! to q, and for which ( 9e) gives 



2y v 8 \b ) E ) 2 


2tt 

IE 


{ 


[A ( 1 + v) - B ( 1 — v )] Sinh 


rr 


b 

l 


+ ( 1 + 


v) 


B rr y Cosh rr 



2ttx 

l 


On posting herein the values (lib) and (lie) for A and B 
it is found that in fact the second summand variable with x 
does not disappear, hence that v( x , ±b/2) does not become 
constant. In order to obtain a value for n independent _of 
k with the least possible deviations from (6b 1 ) 

is put equal to the average value of v at the longitudinal 
plate edge. 

The equation (3) for w can no longer be satisfied, in 
consequence of Rita's formula and is therefore replaced by 
the demand to so choose the free value f that the form 
change energy stored in the plate 

+ 1/2 + b/2 

X= // 

- 1/2 - b/2 

*g *n 

+ 24 ( l -V s ) K A w ) 2 ~ 2 ( 1 - v)(vxx Wyy - v*y)]j-dx dy (12) 
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■becomes a minimum. The unknown f is determined from 

&JL = O’, Differentiation of (12) with respect to f "before 
3 f 

integrating ( Galerkin method), while bearing in mind (la) in- 
stead of (3) gives 


dA 

3f = ~ s 


+ Z/2 + b/2 



$ yy w xx + $ xx w yy ~ 3 ^ xy v xy~ lsf^-v^) AA V 6f dX dy 


1/2 -b/2 


+ b/2 
+ s 

-b/2 


* An + v + w IX s ) 

L yy ^ 5 f Wx af y *y vaf y af/ J 


dy 


x = ± l/2 


+ Z/2 
-Z/2 


du Sw 

a f + Wx a ? 


')] 


d x 


y= ± b/2 


+b/2 

,— jls i_fr' 

12 (1 - V s ) U 
-b/2 

+ Z/2 

*/ 

- 1/2 

+b/2 

+ / | < w xx + u w yy ^ 
-b/2 


( V w xyy - 2 w x yy — w xxx ) 


dy 


aw 
af 

x — ± Z/2 


w xxy “* ^ w xxy ~ w vvv) dW 


dx 


xxy w yyy' 

y = b/2 

dy +j 

x= ± Z/2 -Z/2 


dwx 

af 


( Wyy + VW XX ) 


d x 
y = ± b/2 


+ 2(1- v) 


w xy 


aw 


lb ° 


af 

x = ± z/2 
y = ± b/2 


(13) 
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In equation (13) the "boundary terms and integrals disappear 
with formula (7) and the boundary conditions (6c) and (6d). 

With the solution of f from (13) the stress and strain 
condition of the section can be calculated from equations ( 7) 
and ( 9 ) . 


The subsequent investigation is carried out for the case 
of square bulges (b = l ) , as they occur in the undisturbed 
zone on the sufficiently long cylinder. Entering the ex- 
pression (8) for $ and (11) for 
yields 

e - f 4^ -- A— — ^(2-k} = 0 (14) 


for 

A and B , 

afl 

Ct) (s - K) 

8 


where p is replaced by the expression 

p = s ( £ -t(t7) (15) 

folio w in g from uf±-~-,yj=:be — (cf. equation (9d)) 

and ^ ' 


k = — gg. ? h .... Ln -rJL-- = 0.30935 


tt ( 2 n + S inh 2 tt) 


Posting f = 0 in (14) gives the critical value for the 
shortening e : 




e * 


hence (14) gives 

- 



?( 4 )* - fc (*-«•) ( k - 

2~K 

= 0.5917^) 

( 16) 

The connection between" mean stress 
from ( 15 ) with 

P 

and mean € 

follows 

p - p* = E re a. 

(< - 

€*) 

(17) 
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where 

• Ej. e a = (1 - k) E = i-rJL e = 0.4083 E 

2. - K 

The s tr es s— s tra in equation (17) is reproduced in figure 2. 

The apparent stiffness above the buckling limit is thus only 
about 0.4 times as great as below the critical buckling load. 

, \ 

Entering (16) and (17) in (7) and (9) gives 


2ny 


P* 1 

2 1 


Z* = -JL + —1 ( JL - i) f 

•* l-K \P* ) \ ~ l 

'■ 3 * £ ( a + 2 p) Cosh ^ + 0 SHE sinh^Jj 

ZZ . (JL. - A c 

P* l-K \p* / 


— COS 


2ttx 


1 + a Cosh -~-2- + p —52. Sinh — 


_P_ 
K V p* 


1 sin 


2ttx 


j\a + 0) Sinh -~£ + 0 Cosh 
= — 0 . 


35008, 0 = +0.08429., 


l-K 


= 1.44896 


) 


The numerical evaluation of these equations for simple exceed- 
ing of the buckling load (p/p* - 2) "is shown in figure 3. 


The stress distribution gives only an approximate picture 
of the actual condition, because the deformation of v at the 
section edge is not constant and independent of x as (6b 1 ) 
stipulates. In spite of that the stress distribution is quite 
plain: the wavelike deflections in the side walls after the 
critical buckling load is exceeded tend -to pull the corners 
into the area of maximum deflections. This results in tensile 
stresses in y direction in the middle of the buckles which 
balance the compressive stresses occurring in the reversal 
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point of the waves. This balance can be established only 

with the aid of shearing stresses. The compressive stresses 

in x direction are unevenly distributed after the critical 
buckling load is exceeded, the compressive stresses move 
toward the corners. Depending upon the position of the cross 
section the shearing stresses produce a different type of 
stress distribution over*the plate width so that the shift of 
the compressive stress on the plate edges is greater in the! 
cross section with maximum deflection than in the cross sec- 
tion at the reversal point of the waves. 

The reason that the investigation carried out in the 
foregoing is only approximately valid lies in the form of 
buckle established with formula (7) for which the corners 
were assumed to be square. With increasing exceeding of the 
critical buckling load the corners of the section cease to 
remain square, but pull into the cross sections in which the 
side walls are subjected to the greatest deformations. An 
extension of formula (7) by addition of a second free value 
of the type 

W = T1 + f j cos COS + f s c o 8 cos 

and further terms of this series should certainly produce 
better results. If it is desired to use only one free value 
and to achieve an improvement, formula (7) might be retained 
and other solutions of the homogeneous equation AA $ = 0 
with new constants of integration added in (8b) and these 
then determined by po int— by— po int compliance with the boundary 
conditions. However, the paper work involved is so great that 
its application was foregone. 


EXPERIMENT 


The purpose of the tests was to check the aforementioned 
theoretical results by actual experiments, with special ref- 
erence to formula (7) and the stress— strain formula (17). 

The test specimen was a square cylinder (fig. 4) of 
medium hard sheet brass with the rivet seam .running along the 
center of a side. The characteristics of the brass sheet 
(airplane material identification 2160.7) were: 
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Thickness of sheet s = 0.61 to 0.62 mm, average: 0.615 mm 

Modulus of elasticity E = 10,000 to 11,000 kg/rao 3 , average: 

10,500 kg/mm® 

Tield point a 0>s = 24.4 kg/mm 2 

Tensile strength = 35.2 to 38.6 kg/mm 2 

Elongation 5 = 50 percent 


The wrinkles were measured on two side walls hy means 
of a dial gage and a machined angle section, with arm set 
parallel to two sides of the cylinder as shown in figure 4. 

The sides parallel to the arms were explored hy shifting the 
dial gage. The deformations normal to the direction of com- 
pression were measured in the center of the two side walls 
over the entire length and at two cross sections (test sta- 
tion 18 and 25). The length change of the diagonals under 
increasing load was recorded hy vernier calipers, the compres- 
sion of the test specimen Dy dial gage. The first reading 
was made at 30— kilogram load; the other load stages at which 
measurements were made, were 200, 600, 700, 1000, and 1200 
kilograms. 

The results of the side wall measurements are shown in 
figure 5. According to figures 5a and 5h , the waves formed in 
longitudinal direction are approximately alike so far as the 
uniformity was not disturbed by pre— buckling and edge effects. 
An increase in load was accompanied by a slight change in the 
wave length. As one side area bulged out the area next to it 
bulged in and vice versa. 

The length change of the diagonals is represented in fig- 
ure 6. The absolute size of these deformations is small com- 
pared to the side wall bulging. Serviceable test values start 
at 1300— kilogram loading, shortly before the breaking load is 
reached. The points at which the waves in the side areas have 
their extreme values disclosed a reduction in length of the 
diagonals t , the junction linos of the waves a slight increase 
in corner distance. ' ' 

The result of the compression measurement on the test 
specimen under increasing load is given in table 1, 



HACA TM Ho. 107? 


13 * 


The specimen fails at P = 1330 kilograms through defor- 
mation of the corners at test stations 33 and 34, where the 
waves of the side areas have their junction lines. The failure 
(fig, 7) appears to be instability failure. 

On the basis of the side wall measurements formula (?) 
can be regarded as depicting the deformation area so long as 
the critical load is not too greatly exceeded; if further ex- 
ceeded the corners will not remain square, so that formula 
(7) yields inferior approximate solutions, as already indicated 

With a view to comparing the experimental data for the com- 
pression with the theoretically obtained stress— strain equa- 
tion the values of table 1 were plotted in figure 2, where 

e* = £ _1 (SLBWY = 0.000213 

3 1 - 0 . 3 2 V 8.0 J 

p* = 0.000213 X 1050,000 = 224 kg/cm 2 


Good agreement 
buckling load 
ceeding of the 
elastic range, 
ling load, the 
straight line, 


obtains up to a simple exceeding of the critical 
(p/p* = 2). Failure occurred at 2.3 times ex- 
critical buckling load (p/p* = 3.3). In the 
before reaching the theoretical critical "buck- 
test points do not lie on the theoretical 
which is probably attributable to the pre- 
buckling inevitable on thin sheet. On the basis of the good 
agreement of the theoretical s tr es s—s train curve above the 
critical buckling load the calculation of the constants of 
integration A and B from (lOd) appears therefore justi- 
fied also. The use of (10c) as basis of the study would have 
resulted in a decrease of the apparent stiffness above the 
buckling limit by 0.5 times of the value valid in the elastic 
range, a result obviously less in agreement with the test data 
of figure 2 than that secured from (lOd) with Bred = 0.408 E 
( equation 17 ) , 


A check test of the stress— strain curve was run on a 
specimen of duralumin with the same cross section as that of 
figure 4 (aircraft material 3116.5), (length; J = 800 m/m). 
Its general characteristics were as follows ; 


Sheet thickness s = 0.61 to 0.67 mm,, average = 0,64 no 

Modulus of elasticity E = 6540 to 7320 kg/mm 2 , average 

6830 kg/mn 2 
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Yield limit ff Q>2 = 29.1 kg/an 2 

Tensile strength a B = 41.6 to 42.2 kg/mm 2 


Blongation 6 = 20 percent 


The compression with increasing load 
results are given in table 2, On putting 
table 2 in ratio to the critical values 

«. .Si i , ( P.--M4 V = o 

31-0.3 V 8.0 7 


was measured; the 
the test values of 

000232 


p* ~ 0.000232 X 683,000 # 158 kg/cm 2 


followed by plotting in figure 2 it is found that the theo- 
retical stress— strain curve is here also in good agreement 
with the test points up to a simple exceeding of the critical 
buckling load. Failure occurred at P = 1350 kilograms, 
which corresponds to a 3.1 times exceeding of the critical 
buckling load. 

Summing up it may be stated that the tests fulfilled 
their purpose, namely, to check the formula (?) and- the stress- 
strain equation (17). As a consequence the s tress— strain 
equation (17) secured theoretically with formula (7) can be 
regarded as practical up to a simple exceeding of the critical 
buckling load. 


CONCLUDING NOTE 


On the basis of the present report it is possible to 
obtain a picture of the stress distribution and the strain 
condition in a thin— walled, centrally loaded column of square 
cross section above the critical buckling load, as well as 
an insight into the relationship existing between stress and 
strain. Above the critical buckling load the stress-strain 
curve is linear according to formula (7) in first approxima- 
tion for the square section just as for the stiffened plate 
in compression, except that the range in which this equation 
is applicable, is much great'er compared to the breaking load 
than on the compressed plate. While the test specimens de- 
scribed previously failed at 2.3 and 3.1 times exceeding of 
the critical buckling load, the excess on the stiffened plate 
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in compression may rise to several hundred times above the 
critical buckling load, before the plate fails. The knowledge 
of the str es s— strain equation is important when these sections 
are used as stiffeners on compress ed plates . The decrease in 
section stiffness above the critical buckling load with the 
effective width must then be allowed for. 

The theoretical prediction of the breaking value of thin— 
walled flanged sections in compression is beset with great 
difficulties since, apart from the obstacles indicated in the 
present report, difficulties arise even in the elastic range, 
the stress peaks in the corners exceed the elasticity limit 
before the breaking value is reached and the premises of the 
theory then do no longer hold. If, on the other hand, the 
obtainment of the elastic limit by the stress peaks is re- 
garded as a measure for the load capacity of the sections, the 
methods a°pplied here make it possible to obtain approximate 
values for the allowable loads by theoretical means. 


Translation by J. Vanier, 
National Advisory Committee 
for Aeronautics. 
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TABLE I 


p 

(kg) 

p = 

* r 

(kg /cm 2 ) 

A l 
(mm) 

€ = ' 
l 

0 

0 

0 

0 

200 

102 

0.05 

Ol 

00 

X 

H* 

O 

1 

300 

152 

. 09 

1.05 

400 

203 

.14 

1.63 

500 

254 

.22 

2.56 

600 

305 

.34 

3.95 

700 

355 

.44 

5.11 

800 

406 

.55 

6.39 

900 

457 

.67 

7.79 

1000 

508 

.80 

9.30 

1100 

558 

.92 

10.69 

1200 

6 09 

1.16 

13.48 

1300 

660 

1.47 

17. 08 


TABLE 2 


P 

(kg) 

P 

p “ P 
(kg/ cm 2 ) 

A l 
(mm) 

*-¥ 

' 0 

0 

0 

0 

100' 

49 

.053 

.66 X 10~ 4 

200 

97 

.108 

1.35 •• 

300 

146 

.175 

2.19 

400 

195 

. 29 8 

3.73 

500 

244 

.440 

5,50 

600 

29 3 

.590 

7.38 

700 

341 

.780 

9.75 

800 

390 

.968 

12.10 

900 

439 

1.195 

14,95 

1000 

488 

1.433 

17.90 ; 




Figs. 1,2 
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Figure l.~ Square section under uniform compression. 
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o Test with brass sheet (aircraft material 2160.7) 

» • ii ii duralumin " ( " " 3116.5) 

Theoretical curve according to equation (17) 
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Stresfs strain curves 


Figure 2.- p* = -^ — ("f)^ critical buckling load; e 
critical shortening. 
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Figs. 3,4 
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figures 3a to 3c.- Stress 

distribution 
in the side wall of a thin-wall 



Dial gage holder 


=d centrally loaded section with 
-.quare cross section above the 
iritical buckling load^P^ 



Figure 4.- Dimensions of test 

specimen brass sheet, 
(aircraft material 2160.7); 
dimensions in mm. 
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(5a) Test station I 
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Figures 5a to 5d. Measurements on side walls of test specimen 

(aircraft material 2160.7) 
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Figs. 6,7,8 
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Figure 6.- Length changes in the diagonals (test stations III 
and IV) of the brass 'test specimen (aircraft 
material 2160.7). 


Figure 7.- Specimen after test, 
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Figures 8a and 8b.- Specimen of duralumin (aircraft material 

3116.5) under compressive load of 1200 kg 
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